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Tunable  Oscillations  in  Optically  Injected 
Semiconductor  Lasers  With  Reduced 
Sensitivity  to  Perturbations 

Thomas  B.  Simpson,  Jia-Ming  Liu,  Fellow ;  IEEE ,  Mohammad  AlMulla, 

Nicholas  G.  Usechak,  Senior  Member,  IEEE ,  and  Vassilios  Kovanis 


Abstract — Narrow  line  width  optical  injection  into  a  semiconduc¬ 
tor  laser  can  induce  periodic  oscillations  in  the  injected  laser’s  out¬ 
put  power  with  a  frequency  that  is  widely  tunable  by  simply  varying 
the  steady-state  bias  current  and  operating  temperature.  Recently, 
it  has  been  demonstrated  that  this  oscillation  frequency  can  be 
made  nearly  insensitive  to  small-signal  fluctuations  of  these  two 
parameters  at  certain  operating  points  [1].  Here,  we  demonstrate 
that  this  insensitivity  arises  from  multiwave  mixing  and  interfer¬ 
ence  that  minimizes  the  response  of  both  the  gain  medium  and  the 
circulating  optical  power  at  the  oscillation  frequency.  Both  exper¬ 
imental  measurements  and  model  calculations  of  optical  spectra 
show  that  at  the  operating  points  of  reduced  oscillation  frequency 
sensitivity,  all  strong  components  of  the  optical  spectrum  still  ex¬ 
hibit  a  response  to  the  perturbations.  However,  in  the  power  spec¬ 
tra  and  the  (calculated)  carrier-density  spectra,  the  response  is 
strongly  attenuated.  Novel  operating  points  that  limit  the  sensitiv¬ 
ity  of  the  laser  power  oscillation  frequency  to  perturbations  offer 
the  promise  for  improved  operation  of  tunable  photonic  oscillators 
for  radio-  and  microwave-frequency  applications. 

Index  Terms — Injection  locking,  RF  Oscillators,  tunable  oscilla¬ 
tors. 


I.  Introduction 

SEMICONDUCTOR  lasers,  due  to  their  small  size  and  fast 
dynamical  characteristics,  are  of  interest  not  only  for  their 
use  in  myriad  applications,  but  also  as  a  test  medium  to  un¬ 
derstand  fundamental  characteristics  of  laser  systems  distinct 
from  the  corrupting  influence  of  temperature,  mechanical  and 
other  extrinsic  fluctuations  [2],  [3].  The  injection-locking  and 
dynamical-response  characteristics  of  semiconductor  lasers  sub- 
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ject  to  injected  optical  signals  has  received  considerable  atten¬ 
tion  as  a  tool  for  the  characterization  of  the  intrinsic  dynamic 
parameters  of  semiconductor  lasers  [4]— [8],  and  as  a  mechanism 
for  enhancing  the  fundamental  modulation  characteristics  [9]- 
[12],  with  applications  in  communications  [13] — [16],  RADAR 
[17],  LIDAR  [18],  cryptography  [19],  random  bit  generation 
[20],  and  photonic  frequency  synthesis  [21]-[24].  Under  exter¬ 
nal  injection,  the  semiconductor  laser  can  be  induced  to  display 
a  broad  range  of  output  characteristics  including  stable  injection 
locking,  periodic  power  oscillations  at  the  fundamental  dynamic 
resonance  frequency  of  the  system  and  its  subharmonics,  and 
chaotic  dynamics  [25]-[27].  The  specific  output  characteristics 
can  be  controlled  simply  by  tuning  the  operating  points  of  the 
laser  through  control  of  the  operating  temperature  and  dc-bias 
current.  As  these  capabilities  are  understood,  attention  is  shifting 
to  quantifying  and  minimizing  the  effects  of  unwanted  fluctua¬ 
tions  in  the  operating  point  of  the  laser  (the  dc  bias  current  and 
operating  temperature)  to  the  desired  output  characteristics  [1], 
[23],  [24]. 

A  semiconductor  laser  under  dc  bias  can  have  output  power 
fluctuations  reduced  below  the  standard  quantum  limit  [28]. 
This  has  been  demonstrated  with  an  optically-injected  laser  un¬ 
der  stable  locking  [29].  Less  well  understood  are  the  response 
characteristics  of  optically-injected  lasers  operating  outside  of 
the  stable  locking  regime.  Many  of  the  proposed  applications 
of  the  optically-injected  laser  make  use  of  the  fact  that  tunable 
periodic  or  aperiodic  oscillations  of  the  output  power  can  be 
generated.  Recently,  it  was  demonstrated  that,  in  the  regime 
of  limit-cycle  periodic  dynamics  [the  so-called  period-one  (PI) 
regime],  there  are  operating  points  where  the  PI  oscillation  fre¬ 
quency  could  be  made  essentially  insensitive  to  temperature 
and  bias  current  fluctuations  of  the  injected  laser  [1].  At  specific 
combinations  of  the  bias  point  of  the  injected  laser  and  the  in¬ 
jecting  field’s  amplitude  and  offset  frequency,  the  PI  frequency 
becomes  locally  insensitive  to  the  offset  frequency.  Simulta¬ 
neously,  the  low-frequency  modulation  of  the  slave-laser  bias 
current  is  no  longer  upshifted  to  influence  the  PI  frequency. 
This  is  strictly  a  change  in  the  nonlinear  response;  the  linear 
response  is  largely  unchanged.  It  represents  a  fundamental  shift 
in  small-  signal  sensitivity  of  the  PI  frequency  from  being  dom¬ 
inated  by  changes  in  the  frequency  offset  between  master  and 
slave  lasers  to  being  dominated  by  the  relative  amplitude  of  the 
injected  optical  signal. 

The  objective  of  this  paper  is  to  provide  a  detailed  analy¬ 
sis  of  the  optical  spectra  of  the  optically  injected  semiconductor 


0733-8724  ©  2014  IEEE.  Translations  and  content  mining  are  permitted  for  academic  research  only.  Personal  use  is  also  permitted,  but  republication/redistribution 
requires  IEEE  permission.  See  http://www.ieee.org/publications_standards/publications/rights/index.html  for  more  information. 

1 

Approved  for  public  release;  distribution  unlimited. 


3750 


JOURNAL  OF  LIGHTWAVE  TECHNOLOGY,  VOL.  32,  NO.  20,  OCTOBER  15,  2014 


Fig.  1 .  Schematic  of  the  experimental  apparatus.  VOA,  variable  optical  atten¬ 
uator;  PC,  polarization  controller;  PD,  photodiode.  As  depicted  in  the  figure, 
the  master  and  LO  lasers  were  packaged  with  optical  isolators  while  the  slave 
laser  had  no  isolator.  The  figure  also  identifies  the  propagation  direction  of  the 
light  within  the  setup  using  dashed  arrows. 


laser  undergoing  PI  oscillations.  The  optical  spectra  consist  of  a 
discrete  set  of  strong  spectral  components  and,  it  is  the  complex 
summation  of  the  mixing  of  these  spectral  components  that 
yields  the  optical  power  spectrum.  We  emphasize  this  to  show 
that  at  certain  operating  points,  this  wave  mixing  of  the  mul¬ 
tiple  optical  frequency  components  leads  to  interference  and  a 
cancellation  of  the  response  of  the  PI  oscillation  in  the  power 
spectrum  to  low-frequency  perturbations.  Moreover,  we  quan¬ 
tify  where  this  occurs  in  comparison  to  key  bifurcations  in  the 
dynamics  of  the  laser  system  and  follow  the  changes  of  the  PI 
frequency  at  these  operating  points. 

This  paper  is  organized  as  follows:  Section  II  describes  the 
experimental  configuration  used  to  observe  and  characterize  the 
optical  field  and  power  spectra  of  the  injected  laser.  Section  III 
provides  a  review  of  the  lumped-element  coupled-equation 
model  of  the  optically  injected  semiconductor  laser.  Section  IV 
compares  the  experimental  results  to  those  obtained  from  the 
model,  focusing  on  the  characteristics  of  the  optical  spectra  and 
the  behavior  of  the  fundamental  (PI)  resonance  frequency  of  the 
system.  In  this  section  excellent  quantitative  agreement  between 
this  model  (whose  parameters  have  been  rigorously  determined 
for  this  laser  [26])  and  the  experimental  data  are  obtained.  We 
conclude  the  paper,  in  Section  V,  by  discussing  the  implications 
of  these  results  as  they  pertain  to  applications  for  these  devices. 

II.  Experimental  Apparatus  and  Technique 

The  experimental  apparatus,  shown  schematically  in  Fig.  1, 
is  similar  to  that  previously  used  to  map  the  dynamics  of  a  dis¬ 
tributed  feedback  (DFB)  semiconductor  laser  housed  in  a  butter¬ 
fly  package  containing  an  integrated  thermoelectric  cooler  and 
pigtailed  to  a  single-mode  fiber  but  lacking  an  optical  isolator 
[26] .  By  using  the  same  slave  laser  in  this  paper  we  immediately 
benefit  from  the  experience  obtained  through  previous  studies 
including  apriori  knowledge  of  the  laser’s  parameters.  In  this  pa¬ 
per  the  slave  laser  is  operated  at  three  times  its  threshold  current 
where  it  emits  at  a  wavelength  of  ~1557  nm.  Fight  is  injected 
into  this  laser  from  a  second  (master)  DFB  laser,  this  master 
laser  is  packaged  with  an  optical  isolator  which  suppresses  the 


Fig.  2.  Power  spectrum  of  the  amplified  photodiode  signal  monitoring  the 
output  of  the  optically-injected  laser  operating  in  the  regime  where  its  output 
power  oscillates  at  the  fundamental  resonance  frequency.  The  spectrum  displays 
a  strong  frequency  component  at  this  frequency,  /o ,  and  a  weaker  one  at  its 
second  harmonic,  2/o .  In  this  case  the  normalized  injection  current  was  J  ft*  1.0, 
the  normalized  injection  strength,  £  «  0.02,  and  the  detuning  of  the  master  laser 
was  «  -0.3  GHz. 


same  external  feedback  effects  under  study  in  the  slave  laser. 
A  fiber-based  variable  attenuator  and  an  optical  circulator  are 
placed  between  the  master  and  slave  lasers,  the  former  to  allow 
control  of  the  injected  optical  intensity  and  the  latter  to  increase 
the  optical  isolation  of  the  master  laser  to  >  60  dB .  A  fiber-based 
polarization  controller  is  also  used  to  align  the  polarization  of 
the  master  laser  to  that  of  the  slave.  Both  lasers  are  temperature 
controlled,  allowing  temperature  tuning  of  the  optical  frequency, 
and  both  are  biased  using  precision  current  sources.  As  depicted 
in  the  figure,  weak  radio-  and  microwave-frequency  modulation 
currents  can  be  added  to  the  bias  of  either  the  master  or  slave 
laser  through  bias  Tees.  For  the  experiments  reported  here,  the 
modulation  is  always  weak  (~0.5%  peak  amplitude)  relative 
to  the  dc  bias,  and  is  fixed  at  a  frequency  of  500  MHz.  The 
output  of  the  slave  laser  is  detected  by  a  high-speed  photodi¬ 
ode  whose  output  is  subsequently  amplified  and  then  monitored 
by  an  electronic  spectrum  analyzer  (ESA).  With  just  the  slave 
laser  output,  the  power  spectrum  of  the  laser  output  is  gener¬ 
ated.  Alternatively,  the  slave  laser  output  can  be  combined  in  a 
fiber  coupler  with  the  output  from  a  third,  temperature-tunable 
DFB  laser  that  acts  as  a  local  oscillator  (FO).  In  this  heterodyne 
case,  the  electrical  spectrum  of  the  photodiode’s  output  serves 
as  a  high-resolution  measurement  of  the  amplitude  spectrum  of 
the  optical  field.  Use  of  the  heterodyne  spectrum  yields  a  much 
higher  resolution,  limited  by  the  ~1  MHz  linewidth  of  the  FO 
laser,  than  can  be  achieved  using  either  a  conventional  optical 
spectrum  analyzer  or  a  high-finesse  Fabry-Perot  cavity. 

Under  optical  injection,  the  slave  laser  output  can  become 
destabilized,  which  can  in  turn  result  in  power  oscillations.  Fig.  2 
shows  an  example  of  the  power  spectrum  generated  by  the  laser 
under  optical  injection  in  the  PI  oscillation  regime.  It  is  domi¬ 
nated  by  a  strong  peak  at  the  fundamental  resonance  frequency 
(/o ,  ^6  GHz)  and  weaker  peaks  at  the  harmonics  of  this  fre¬ 
quency.  When  the  laser  generates  an  output  with  periodic  power 
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Fig.  3.  Heterodyne  measurement  of  the  optical  frequency  spectra  generated 
by  sweeping  the  LO  laser  frequency  across  the  injected  laser  output  spectrum 
by  ramping  the  operating  temperature.  The  offset  frequency  is  with  respect  to 
the  free-running  laser  output  frequency,  vco .  Peaks  are  labeled  by  their  optical 
frequency.  Shown  are  the  spectrum  of  the  free-running  laser  (red)  with  its 
weak,  but  broad,  relaxation  resonance  frequency  sidebands,  and  the  spectrum 
of  the  optically  injected  semiconductor  laser  operating  in  the  region  of  power 
oscillations  at  the  fundamental  resonance  frequency.  In  the  latter  case,  the  optical 
injection  frequency,  vY,  was  at  the  frequency  corresponding  to  the  peak  just 
below  the  optical  frequency  of  the  free-running  laser,  and  the  pushed  laser 
cavity  resonance  frequency  is  vc  =  vx  —  /o .  The  operating  conditions  here  are 
the  same  as  in  Fig.  2. 


Fig.  4.  A  comparison  of  two  ways  of  generating  the  optical  frequency  spec¬ 
trum  through  heterodyne  mixing  with  a  local  oscillator.  The  left  panel  shows 
two  spectra  generated  by  sweeping  the  LO  laser  through  the  spectrum  of  the 
optically  injected  laser  as  in  Fig.  3.  The  right  panel  shows  two  spectra  generated 
when  the  LO  laser  frequency,  z^lo  ,  is  fixed  at  the  optical  frequency  indicated  by 
the  green  arrow  in  the  left  panel.  The  green  trace  shows  the  spectrum  generated 
when  the  slave  laser  bias  current  is  turned  off,  showing  the  beat  signal  gener¬ 
ated  by  the  mixing  of  the  LO  laser  with  the  residual  master  laser  output  that  is 
scattered  off  of  the  front  facet  of  the  slave  laser  and  gives  the  relative  frequency 
of  the  master  laser.  The  black  trace  shows  the  dominant  peaks  of  the  combined 
power  and  heterodyne  spectrum,  with  the  optical  spectrum  peaks  positioned 
at  their  relative  offsets  with  respect  to  the  LO  laser  frequency.  For  the  injected 
laser  there  is  a  strong  peak  at  the  injection  frequency  {yx),  at  the  injection- shifted 
laser  cavity  oscillation  frequency  ( vc )  that  is  offset  by  the  PI  frequency  deter¬ 
mined  by  the  strong  peak  in  the  power  spectrum  (/o),  and  weaker  peaks  offset 
by  further  PI  frequency  shifts.  The  free-running  laser  spectrum  peak  at  z/co  is 
shown  in  the  left  panel  but  not  the  right  panel.  Here  the  master-laser  offset  was 
—4.5  GHz,  a  normalized  injection  strength  of  £  =  0.06  was  used,  and  the  slave 
laser  was  biased  three-times  above  its  threshold  current. 


oscillations  under  optical  injection,  it  emits  at  a  discrete  set 
of  optical  frequencies  separated  by  the  characteristic  resonance 
frequency.  Fig.  3  shows  the  optical  spectrum  of  the  injected  laser 
in  comparison  to  that  of  the  free-running  laser.  The  spectra  are 
generated  by  fixing  the  measurement  frequency  of  the  ESA  and 
slowly  ramping  the  temperature,  and  thus  the  wavelength,  of  the 
LO  laser.  When  the  beat  frequency  of  the  LO  and  one  of  the 
field  components  of  the  slave  laser  matches  the  fixed  frequency 
of  the  ESA,  a  strong  signal  is  detected.  Typically,  there  are  2  to 
4  dominant  frequency  peaks  of  the  slave  laser  when  operated  in 
the  PI  oscillation  regime  with  the  two  strongest  spectral  compo¬ 
nents  at  the  injection  frequency,  v\,  and  at  the  injection- shifted 
laser-cavity  resonance  frequency,  isc  =  —  /o .  For  clarifica¬ 

tion  purposes  we  point  out  that,  to  improve  readability,  this 
paper  uses  v  to  define  optical  frequencies  whereas  /  is  used 
to  define  microwave  frequencies.  The  free-running  laser  has  a 
single  dominant  peak  at  the  laser  cavity  resonance  frequency, 
vc o,  and  weak  broad  sidebands  peaking  at  approximately  the 
relaxation  resonance  frequency.  The  ESA  is  set  to  a  value  of 
50  MHz  to  generate  the  spectra  of  Fig.  3.  This  yields  relatively 
low-resolution  optical  spectra,  due  to  the  double- valued  nature 
of  the  beat  frequency  (either  positive  or  negative  offset  generates 
a  signal).  Although  they  obscure  detailed  characteristics  of  the 
optical  frequency  components,  such  as  changes  due  to  radio¬ 
frequency  current  modulations,  these  spectra  are  very  useful  for 
mapping  out  the  dynamic  operating  characteristics  of  the  laser. 

To  achieve  higher  resolutions,  the  optical  frequency  of  the 
LO  laser  can  be  fixed.  In  this  study,  we  are  specifically  in¬ 
terested  in  the  characteristics  of  the  optically  injected  laser 
when  its  output  is  undergoing  periodic  power  oscillations  with  a 


tunable  period-one  (PI)  resonance  frequency.  In  this  case  the 
strong  spectral  components  of  the  PI  oscillation  are  separated 
by  several  gigahertz;  therefore,  the  LO  laser  can  be  tuned  to  an 
optical  frequency  between  the  spectral  peaks.  Here,  the  spec¬ 
trum  of  the  photodiode  output  consists  of  the  combined  power 
spectrum  of  the  slave  laser  mixing  with  itself  and  the  hetero¬ 
dyne  spectrum  of  the  mixing  of  the  two  lasers.  Because  the  LO 
laser  output  is  essentially  monochromatic,  its  power  spectrum 
has  no  spectral  features  that  contribute  to  our  optical  spectrum 
or  degrade  our  results.  An  example  demonstrating  this  claim  is 
shown  in  Fig.  4.  In  this  figure  the  left  panel  displays  the  spectra 
generated  by  tuning  the  LO  laser  (as  done  in  Fig.  3)  and  iden¬ 
tifies  the  frequency  where  the  LO  laser  was  fixed  to  generate 
the  spectra  shown  in  the  right  panel.  While  care  must  be  taken 
in  the  positioning  of  the  LO  frequency,  it  is  possible  to  inter¬ 
pret  the  spectra  to  identify  the  optical  frequency  components 
at  positive  and  negative  offsets  with  respect  to  the  LO  and  the 
power  spectrum  peak  at  the  PI  frequency.  The  spectral  position 
of  the  peak  at  the  master  laser  frequency,  is  identified  by 
turning  off  the  slave  laser  and  using  the  small  reflected  signal  of 
the  master  laser  off  the  slave  laser  facet  as  shown  in  the  figure. 
Similarly,  the  spectral  position  of  the  free-running  laser,  uco, 
is  established  by  turning  off  the  master  laser,  and  the  position 
of  the  PI  peak,  /0,  in  the  power  spectrum  is  established  by 
turning  off  the  LO  laser.  With  this  technique  the  resolution  is 
limited  by  the  1-MHz  linewidth  of  our  LO,  provided  the  spectral 
features  do  not  overlap.  Each  of  the  prominent  optical  features 
of  the  slave  laser  can  be  distinguished  by  properly  positioning 
the  output  frequency  of  the  LO  between  the  emission  peaks 
of  the  slave.  In  our  laboratory  setup  both  the  master  and  LO 
lasers  are  of  the  same  commercial  type  and  display  a  nearly 
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Fig.  5.  A  spectrum  similar  to  the  right  panel  of  Fig.  4  but  now  showing  the 
modifications  due  to  adding  a  weak  current  modulation  at  500  MHz  to  the  bias 
of  the  optically  injected  laser.  Peaks  labeled  /m  and  fo,  and  their  associated 
lower  and  upper  sidebands,  are  the  peaks  at  the  modulation  and  PI  frequency, 
respectively,  in  the  power  spectrum.  Peaks  labeled  with  v\ ,  uc,  +  fo,  and 
v i  —  fo,  with  associated  sidebands,  are  from  the  optical  spectrum.  See  text 
for  further  identification  of  all  of  the  peaks  in  the  heterodyne  spectrum.  The 
frequency  axis  is  as  measured  by  the  ESA  and  so  the  positions  of  the  heterodyne 
peaks  from  the  optical  spectrum  depend  on  the  offset  with  respect  to  the  local 
oscillator  laser  frequency,  •  Here  the  master-laser  offset  was  —3.1  GHz, 
a  normalized  injection  strength  of  £  =  0.06  was  used,  and  the  slave  laser  was 
biased  three-times  above  its  threshold  current. 


identical  tuning  characteristic  with  temperature  of  —11.5 
GHz/C,  which  makes  it  relatively  straightforward  to  keep  the 
LO  laser  frequency  properly  positioned  relative  to  the  slave  laser 
output  spectrum  as  the  master  laser’s  optical  frequency  is  tuned. 

In  this  paper,  we  add  a  weak  modulation  current  to  the  dc 
bias  of  either  the  master  or  slave  laser  in  order  to  investigate  the 
resulting  spectral  properties.  This  modulation  adds  negatively 
and  positively  offset  sidebands  to  each  of  the  strong  optical- 
field  components,  as  shown  in  Fig.  5  ,  where  all  of  the  labeled 
peaks  can  be  identified.  From  the  power  spectrum  there  are 
peaks  at  the  500-MHz  modulation  frequency  (labeled  /m)  and 
at  the  PI  frequency  and  its  upper  and  lower  sidebands,  labeled 
/o ,  fo  +  /m ,  and  fo  —  fm  respectively.  Additionally,  the  het¬ 
erodyne  spectrum  adds  peaks  due  to  mixing  between  the  LO 
(u lo  )  and  the  components  and  sidebands  that  are  associated  with 
the  master  laser  injection  frequency  (yx,  vx  +  /m,  and  v\  —  /m), 
the  injection- shifted,  slave-laser  cavity  resonance  frequency  (uc , 
uc  +  fm,  and  z'c  —  /m),  the  component  one  PI  frequency  offset 
above  the  master  laser  ( v\  +  fo),  and  the  component  one  PI  fre¬ 
quency  offset  below  the  pushed  frequency  peak  (z/c  —  fo).  (Note 
that  the  data  in  the  right  panel  of  Fig.  4  and  in  Fig.  5  are  the  raw 
data  from  the  ESA.  The  values  taken  from  such  spectra  and  used 
in  subsequent  figures  have  been  corrected  for  the  response  char¬ 
acteristics  of  the  measurement  system.)  The  amplitudes  of  all 
of  the  sidebands  vary  in  nontrivial  ways  as  a  control  parameter 
such  as  the  frequency  detuning  between  the  master  laser  and  the 
free-running  slave  laser  is  varied.  The  sidebands  of  each  optical 
field  component  mix  with  the  adjacent  central  feature  of  the 
adjacent  field  components  to  produce  the  sidebands  on  either 
side  of  the  PI  frequency  feature  in  the  power  spectrum.  This 


complex  summation  and  mixing  can  lead  to  destructive  inter¬ 
ference,  a  feature  that  can  be  recovered  in  the  standard  lumped- 
element  coupled-equation  model  of  the  injected  laser  system 
that  we  review  next. 


III.  Theoretical  Model 


The  output  from  a  laser  subject  to  an  externally  injected  field 
can  be  described  by  a  set  of  coupled  equations  relating  the  cir¬ 
culating  optical  field,  the  polarization  and  the  population  den¬ 
sity  of  the  gain  medium,  and  the  injected  optical  signals  [30]. 
Specifically,  for  the  case  of  a  semiconductor  laser  subject  to 
near-resonant  external  optical  injection  the  coupled  equations 
for  the  circulating  optical  field,  A{t)e~luJcot ,  where  cjco  is  the 
optical  frequency  of  the  free-running  laser,  and  the  carrier  den¬ 
sity  of  the  gain  medium,  N(t)  simplify  to  [26] 


dA 

dt 


~  ■“  +  i(^c0  —  idcc) 


+  -g{N, \A\)A  +  t)A  ie-i(wi-Wc»)<: 
d,N  J  2  e0n2  ,^T  .  ....  . 

def 


g{NM\)  =  Re[g(N,\A\)\ 
g(N,\A\)  =  g0(N0,A0)+ 


(1) 

(2) 

(3) 


(N-No) 


N0,Ao 


N0,Ao 

(\A\2~A20).  (4) 


Here,  the  semiconductor  laser  is  modeled  as  a  single-mode 
lumped  element  subject  to  an  injected  laser  beam  of  con¬ 
stant  amplitude,  A[,  and  constant  optical  frequency  cji  =  27rz/i. 
The  laser  cavity  is  described  in  terms  of  7C,  the  cavity  decay 
rate,  cuco  =  27 rz/c0 ,  the  free-running  cavity  oscillation  frequency, 
c jcc  =  27rucc,  the  cold-cavity  (carrier-free)  resonance  frequency, 
and  T,  the  confinement  factor  that  relates  the  overlap  of  the  os¬ 
cillating  mode  with  the  gain  medium,  The  injection  coupling 
rate,  g,  quantifies  the  coupling  of  the  externally  injected  laser 
beam  into  the  slave  laser  cavity.  The  carrier  density  injection 
rate  is  J/ed  where  J  is  the  injected  current  density  ( e  is  the 
electronic  charge  and  d  is  the  thickness  of  the  gain  medium),  7S 
is  the  carrier  spontaneous  decay  rate,  eon2  is  the  non-resonant 
(real)  permittivity  of  the  semiconductor  gain  medium,  and  h  is 
Plancks  constant.  When  free-running,  in  other  words  with  no 
injected  laser  beam,  the  laser  operates  at  an  optical  frequency 
cjco  with  a  steady-state  amplitude  Ao  and  a  steady-state  carrier 
density  No.  The  complex  gain,  g(N,  \A\),  controls  the  charac¬ 
teristics  of  the  laser  cavity  and  depends  on  the  carrier  density  and 
the  magnitude  of  the  laser  amplitude.  The  real  part  of  the  gain, 
g(N,  |  A  | ) ,  describes  the  transfer  of  power  between  the  semicon¬ 
ductor  gain  medium  and  the  laser  beam,  while  the  imaginary  part 
describes  how  the  dielectric  constant  (and  through  it  the  laser 
cavity  resonance  frequency)  changes  in  the  laser.  Both  parts 
are  assumed  to  vary  about  the  free-running,  steady-state  oper¬ 
ating  point,  go  (Nq ,  |  A0 1),  as  a  linear  function  of  the  changes  in 
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the  carrier  density  and  the  laser  intensity  (proportional  to  the 
square  of  the  laser  field  magnitude).  In  (4),  b  is  the  ratio  of 
the  imaginary  to  real  parts  of  the  complex  gain  due  to  carrier 
density  changes,  and  b'  is  the  ratio  due  to  saturation  effects. 
The  latter  ratio  is  not  generally  used  in  describing  semiconduc¬ 
tor  lasers,  but  is  essential  for  the  quantitative  modeling  done 
here  [5],  [26].  Just  as  b  was  originally  introduced  to  account 
for  the  enhanced  linewidth  of  semiconductor  lasers,  and  satura¬ 
tion  effects  (related  to  ^w)  were  introduced  to  account  for  the 
low-frequency  rolloff  in  the  modulation  response,  b'  was  intro¬ 
duced  to  account  for  the  asymmetry  of  the  optical  spectrum  in 
response  to  a  modulation  current  [5].  It  is  generally  expected  to 
be  present  in  a  detuned  two-level  system  [3 1] ,  [32] .  Equation  (4) 
gives  the  general,  first-order  expansion  of  the  gain  in  terms  of 
changes  to  the  carrier  density  and  optical  intensity;  it  does  not 
rely  on  a  specific  quantum  model  of  the  gain  medium.  It  as¬ 
sumes  that  the  carriers  involved  in  the  optical  transition  come 
into  equilibrium  with  all  carriers  on  times  scales  short  com¬ 
pared  to  those  of  interest  for  the  dynamics.  This  is  expected  to 
be  true  of  most  semiconductor  laser  media  and  the  model  has 
been  successfully  applied  to  a  wide  range  of  lasers  based  on  dif¬ 
ferent  cavity  configurations  (Fabry-Perot,  DFB,  VCSEF,  etc.) 
and  materials  (heterostructure,  quantum  well,  and  a  variety  of 
III-V  alloys).  However,  recent  work  on  quantum-dot  semicon¬ 
ductor  lasers,  and  indeed  the  whole  rationale  for  their  use,  in¬ 
dicates  that  the  quantum-confined  carriers  involved  in  the  opti¬ 
cal  transition  need  to  be  considered  separately  from  the  back¬ 
ground  carriers  in  the  non-localized  energy  bands  [33],  [34]. 
The  appropriate  lowest-order  model  for  this  gain  medium  is 
still  the  subject  of  debate  and  ongoing  research.  More  com¬ 
plex  laser  structures,  such  as  multi- section  lasers,  could  still 
conceivably  be  treated  as  a  lumped  element  but  all  sections 
would  require  appropriate  equations  for  their  respective  carrier 
dynamics,  and  the  appropriateness  of  the  single- spatial  mode 
approximation  would  need  to  be  determined.  Finally,  we  do  not 
include  spontaneous  emission  in  the  work  here  because  in  our 
experiments  the  external  optical  injection  and  deterministic  in¬ 
teractions  are  strong  enough  to  dominate  the  dynamics  of  the 
problem.  The  spontaneous  emission  plays  a  role  in  determin¬ 
ing  the  linewidth  of  the  various  peaks  in  the  optical  spectrum 
(excepting  the  peak  at  the  injection  frequency),  but  otherwise 
acts  as  a  perturbing  noise  source  that  sets  the  noise  floor  for  the 
spectra. 

Equation  (1)  is  the  equation  for  the  complex  field,  A.  For 
numerical  modeling  it  is  convenient  to  introduce  A{t)e~lu;cot  = 
A0e~luJit  [ai(t)  A  ia2(t)]  and  TV  =  7V0  [1  +  n\  to  generate  two 
real  equations  for  the  field  components  in-phase  and  in¬ 
quadrature  with  respect  to  the  master  laser  ([1]  and  its  sup¬ 
plementary  material) 


da\ 

dt 


—Cla2  +  7c£  + 
7p  /  2  I  2 

—  —  [CLi  +  a2 

7c 


7c  7n  ~  /  .  ,  x 

—  — ~n  (ai  +  ba2) 

2  l_7s J 

—  l)  (ai  +  b'a2) 


(5) 


TABLE  I 

Parameter  Values  Used  in  Numerical  Simulations 


Term 

Description 

Value 

Units 

Error 

b 

Linewidth  Enhance.  Factor 

3.0 

N/A 

+  10/— 20% 

b' 

Nonlinear  Gain  Ratio 

-1.5 

N/A 

±50% 

7c 

Cavity  Photon  Decay  Rate 

5.4  x  1011 

s-1 

±10% 

7s 

Carrier  Spont.  Decay  Rate 

4  x  109 

s-1 

±10% 

7n 

Differential  Gain  Relax.  Rate 

1.7  x  10 9  •  J 

s-1 

±10% 

7P 

Nonlinear  Gain  Relax.  Rate 

5.0  x  10 9  •  J 

s-1 

±10% 

^  =  S)«,  +  7 


dt 


(—bai  A  a2) 


2  HJ 


—  —  (a\  +  a2  —  l)  (—b'ai  A  a2) 

7c 


In  these  equations,  7n  =  gjL 


A0,N0 


(6) 

At  7P  = 


-r 


dg 


d\A[ 


Ao  ,Nq 


Ao,  and  J  =  -  %Nq)  /(7sjV0)  where 

Jo  is  the  steady-state  injection  current,  £  =  77  and 
Cl  =  uj[  —  cjco  .  The  equation  for  the  carrier  density  can  be 
written  as 


dn 

dt 


—  7s  (j  A  l)  -j - [7s  +  7n (al  A  al)] 


—  If sJ 


Jo 

l_>(a2+a2} 

7c 


{a2,  +a22-l)  (7) 


where  Jm  is  any  time-dependent  part  of  the  injection  current. 
The  laser  parameters  can  be  experimentally  determined  [5],  [6] 
while  J,  Cl,  and  £  can  be  varied  in  the  experiment.  The 
model  calculations  use  the  experimentally  determined  parame¬ 
ters  listed  in  Table  I  for  the  injected  laser  [1]  biased  at  J  ~  2.0. 
In  these  experiments  ^-<1. 

For  the  free-running  laser,  the  relaxation  resonance  frequency 
is  given  by  A  7s7P  and  the  associated  damping  rate  by 

7s  +  7n  +  7P  •  Note  that  because  is  approximately  two  or¬ 
ders  of  magnitude  larger  than  the  other  rates,  the  relaxation 
resonance  frequency  is  typically  much  larger  than  the  damping 
rate  and  is  only  weakly  dependent  on  the  nonlinear  gain  relax¬ 
ation  rate.  Because  of  this,  one  can  get  good  semi-quantitative 
agreement  between  model  and  experiment  even  with  7P  =  0, 
or  equivalently  =  0  [35].  As  a  consequence,  all  of  the 
key  experimentally  measured  features  of  the  optically  injected 
semiconductor  laser  discussed  in  this  paper  can  be  reproduced 
with  7P  =  0  (and  a  modified  value  of  b),  but  at  significantly 
offset  values  of  the  detuning  frequency  and  injection  parameter 
relative  to  the  experimentally  measured  values. 

Numerical  simulations  of  the  optically  injected  semiconduc¬ 
tor  laser  based  on  (4)— (7)  used  a  Runge-Kutta  integration  pro¬ 
cedure  over  durations  >  1/rs  to  generate  time  series  that  were 
subsequently  Fourier  transformed  to  produce  spectra  [13]  that 
could  be  directly  compared  with  the  experimentally  measured 
spectra.  Fig.  6  shows  a  calculated  mapping  of  the  dynamics 
regions  overlaid  with  the  curves  of  constant  PI  frequency. 
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Injection  Strength 

Fig.  6.  Mapping  of  the  dynamics  of  the  optically  injected  semiconductor  laser 
as  a  function  of  the  strength  and  offset  frequency  of  the  master  laser.  The  white 
region  is  stable  injection  locking,  the  large  multicolored  region  displays  power 
oscillations  at  the  PI  frequency,  with  lines  of  constant  oscillation  frequency 
displayed.  Dark  blue  to  black  are  regions  of  period  doubling  of  the  fundamental 
resonance  and  more  complex  dynamics.  Also  drawn  is  a  vertical  line  locating 
the  operating  conditions  for  data  displayed  in  subsequent  figures. 


The  region  of  PI  dynamics  as  a  function  of  the  master  laser 
detuning  frequency  is  bounded  from  below  by  a  Hopf  bifurca¬ 
tion  to  stable  locking,  and  in  certain  regions  by  transitions  to 
more  complex  periodic  (period-doubling  and  higher)  and  ape¬ 
riodic  (quasiperiodic  and  chaotic)  dynamics.  This  mapping  is 
in  quantitative  agreement  with  the  experimentally  observed  dy¬ 
namical  characteristics  of  the  laser.  Note  that  there  is  a  line 
of  operating  points  extending  from  near  the  minimum  offset 
of  the  Hopf  bifurcation  to  the  upper  edge  of  lower  amplitude 
and  offset  region  of  complex  dynamics  where  the  PI  frequency 
becomes  locally  insensitive  to  changes  in  the  offset  frequency 
between  master  and  slave  lasers.  We  will  emphasize  these  oper¬ 
ating  points  below,  specifically  detailing  how  the  optical  spectra 
change  as  the  master  laser  detuning  is  varied  in  the  regions 
around  these  operating  points. 

IV.  Results:  Experimental  Measurements 
and  Numerical  Simulations 

As  one  moves  around  the  map  of  Fig.  6,  by  changing  the 
injected  amplitude  (varying  the  master  laser  bias  current  or  the 
attenuation  between  the  master  and  slave  lasers)  or  offset  fre¬ 
quency  (varying  the  operating  temperature  of  the  master  laser), 
changes  in  the  optical  heterodyne  spectra  and  the  power  spec¬ 
tra  are  observed,  even  within  the  region  of  PI  dynamics.  For 
instance,  as  one  moves  along  the  line  in  Fig.  6  away  from  the 
Hopf  bifurcation  separating  the  PI  dynamics  region  from  the  re¬ 
gion  of  stable  locking,  the  strong  peaks  in  the  optical  spectra  and 
the  power  spectra  change  in  amplitude  and  frequency  as  shown 
in  Fig.  7.  The  optical  peak  at  the  injection  frequency  decreases 
in  amplitude,  while  the  other  optical  peaks  first  grow  stronger, 
then  level  off  and  finally  begin  to  decrease  in  amplitude.  The  PI 
frequency  first  decreases,  passes  through  a  minimum  and  then 
increases.  Note  that  the  offset  for  the  PI  minimum  is  close  to 


Offset  Frequency  (GHz) 

— ★ —  Master  frequency,  vj}  peak  • - Above  master,  Vj+f0,  peak 

— ^ —  Pushed  frequency,  vc,  peak  — ▼ - Below  pushed,  vc-f0,  peak 

- ♦ - PI  frequency,  f0,  peak  B - PI  frequency 

Fig.  7.  Experimentally  measured  amplitudes  of  the  various  optical  frequency 
components,  and  the  component  in  the  power  spectrum  at  the  PI  resonance 
frequency,  /o,  as  the  master  laser  detuning  is  varied  along  the  vertical  line 
superimposed  in  Fig.  6.  The  master  frequency  peak  is  at  the  injection  frequency 
of  the  master  laser,  vY ,  the  pushed  laser  cavity  peak  is  offset  by  the  PI  frequency 
below  this  value,  isc  —  v-x  —  /0 ,  and  the  other  two  are  for  the  peaks  above 
(vi  +  /o)  and  below  (uc  —  /o)  these  central  peaks.  Also  shown  is  the  value  of 
the  PI  frequency  (right  axis). 
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Fig.  8.  Mapping  of  key  operating  points  experimentally  measured  in  this 
investigation.  See  text  for  details.  SNLC  =  saddle  node  of  limit  cycle  bifurcation, 
Master  Mod  =  minimum  sensitivity  for  modulation  current  added  to  the  master 
laser.  The  slave  laser  modulation  minima  overlap  with  the  PI  minima  that  run 
from  the  SNLC  bifurcation  to  the  Hopf  bifurcation.  The  normalized  injection 
power  (=  10.0)  is  at  the  same  injection  strength  as  £  =  0.01. 


the  offset  where  the  two  dominant  optical  peaks,  at  the  injection 
frequency  and  at  the  pushed  cavity  resonance  frequency  that  is 
offset  by  the  PI  frequency,  are  of  the  same  amplitude. 

Operating  points  where  the  PI  frequency  is  at  an  extremum 
with  respect  to  changes  in  a  control  parameter  are  of  particular 
interest  because  they  show  local  insensitivity  to  perturbations  in 
that  parameter.  Because  of  the  application  potential  for  these  re¬ 
gions,  we  set  out  to  find  such  extrema  and  better  understand 
them.  These  results  are  summarized  in  Fig.  8  which  shows 
a  mapping  of  local  sensitivity  minima  with  respect  to  some 
variation  along  with  key  bifurcations  as  the  detuning  between 
master  and  slave  laser  is  varied.  Starting  from  negative  offsets 
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Fig.  9.  The  PI  resonance  frequency  just  above  the  Hopf  and  saddle-node 
of  limit  cycle  bifurcations  shown  in  Fig.  8,  and  at  the  operating  points  where 
the  PI  frequency  and  the  sensitivity  to  current  fluctuations  are  simultaneously 
minimized.  The  saddle-node  of  limit  cycles  bifurcation  ends  when  the  minimum 
PI  frequency  intersects.  The  normalized  injection  power  (=  10.0)  is  at  the  same 
injection  strength  as  £  =  0.01. 


Fig.  10.  Experimentally  measured  (symbols)  and  calculated  (solid  curves) 
values  of  the  modulation  sideband  amplitude  and  the  PI  frequency  using  a  nor¬ 
malized  injection  amplitude  of  £  =  0.06  when  the  current  modulation  is  applied 
to  the  slave  laser.  This  figure  shows  the  positive  frequency  detuning  sideband 
(solid  circles,  corresponding  to  peak  /o  +  /m  in  Fig.  5),  the  negative  detuning 
sideband  (solid  triangles,  corresponding  to  peak  /o  —  /m  in  Fig.  5),  and  the 
corresponding  PI  frequency  (open  squares,  corresponding  to  the  frequency  of 
the  peak  labeled  /o  in  Fig.  5). 


and  working  up,  the  region  of  stable  locking  is  bounded  by 
saddle-node  (open  squares;  purple)  and  Hopf  (sold  triangles; 
green)  bifurcations.  Above  the  Hopf  bifurcation  is  the  region  of 
limit-cycle  oscillations  at  the  PI  frequency.  There  is  a  bounded 
range  of  chaotic  dynamics  that  from  below  shows  period¬ 
doubling  routes  to  chaos  (open  circles;  orange),  and  from  above 
shows  an  abrupt  saddle-node  of  limit  cycle  (SNLC)  bifurca¬ 
tion  (solid  circles;  orange)  [25].  On  either  side  of  the  region  of 
chaotic  dynamics  there  are  local  minima  in  the  PI  frequency.  At 
low  injection  levels  a  range  (solid  diamonds;  pink)  runs  down 
to  near  the  upper  frequency  bound  on  the  Hopf  bifurcation.  On 
the  high  injection  side  of  the  chaotic  dynamics  region  there  are 
two  ranges,  a  short  one  associated  with  a  period-doubling  route 
from/to  chaos  (solid  right-facing  triangles;  gray)  and  a  more 
extensive  one  that  runs  between  the  SNLCs  bifurcation  and  the 
Hopf  bifurcation  (solid  diamonds;  brown).  Fig.  9  shows  how 
the  frequencies  of  the  PI  minima  compare  with  those  along 
the  Hopf  and  SNLC  bifurcations  as  a  function  of  the  injection 
strength.  In  Fig.  9,  the  SNLCs  terminates  when  its  frequency 
matches  the  PI  minimum  frequency  at  the  same  optical  injec¬ 
tion  level  [36].  If  the  laser  bias  is  changed,  the  PI  minima  can 
terminate  in  the  region  of  complex  dynamics  close  to  this  bifur¬ 
cation.  In  addition  to  minima  in  the  PI  frequency  with  respect 
to  detuning,  we  searched  for  minima  in  the  sensitivity  of  the  PI 
frequency  to  modulation  of  the  bias  current  to  either  the  master 
or  slave  laser.  For  modulation  of  the  slave  laser  current,  we  ob¬ 
served  that  the  sensitivity  was  also  minimized  along  the  curve 
of  brown  diamonds  in  Fig.  8,  coincident  with  the  PI  minimum 
with  respect  to  detuning.  Fig.  10  shows  the  dramatic  decrease  in 
response  of  the  sidebands  around  the  PI  frequency  peak  as  the 
master- slave  frequency  offset  is  varied  along  the  line  denoted 
in  Fig.  6.  For  modulations  to  the  master  laser  current,  there  is  a 
separate  range  denoted  by  the  solid  blue  squares  (Master  Mod) 
in  Fig.  8  that  are  offset  from  the  PI  frequency  minima.  These 
minima  are  offset  because  changes  to  the  master  laser  current 


change  the  amplitude  of  the  injected  optical  signal  to  the  slave 
laser  directly  modifying  the  circulating  optical  field,  a  different 
perturbation  to  the  laser  system  than  a  current  modulation  to 
the  interaction  gain  medium  that  directly  modifies  the  carrier 
density. 

A  close  inspection  of  Fig.  10  shows  that  the  minima  of  the 
two  sidebands  are  offset  with  respect  to  each  other  and  to  the  PI 
frequency  minimum  with  respect  to  detuning.  This  is  clearest 
in  the  numerical  data,  which  is  not  corrupted  by  noise  or  other 
experimental  artifacts.  The  offset  with  respect  to  each  other  is  a 
monotonically  increasing  function  with  respect  to  the  modula¬ 
tion  frequency.  Calculations  show  that  the  offset  with  respect  to 
the  PI  frequency  minimum  is  a  complex  function  of  the  various 
dynamic  parameters,  in  particular  the  parameter  b'  [31]. 

The  set  of  operating  points  where  there  is  the  simultaneous 
minimum  of  the  PI  frequency  with  respect  to  detuning  and 
sensitivity  minimum  to  current  modulation  of  the  slave  laser  is 
of  particular  interest  due  to  the  fact  that  it  represents  operation 
where  sensitivity  to  current  and  temperature  fluctuations  of  the 
slave  laser  are  minimized,  and  simultaneously  where  sensitivity 
to  temperature  (but  not  current)  fluctuations  of  the  master  are 
minimized.  Therefore,  we  investigated  the  changes  to  the  optical 
spectra  along  the  detuning  line  of  Fig.  6  in  greater  detail  in 
Fig.  11.  Here,  the  sideband  signal  amplitude,  normalized  to 
the  amplitude  of  its  associated  center  peak,  is  shown  as  the 
master  laser  detuning  frequency  is  varied.  The  sidebands  of 
each  of  the  four  strongest  peaks  in  the  optical  spectrum  display 
different  characteristics  as  the  detuning  is  changed.  In  no  case 
is  the  minimum  among  the  sidebands  of  the  various  optical 
peaks  as  deep  as  the  sensitivity  minimum  of  the  sidebands  in 
the  power  spectrum.  Moreover,  the  linear  modulation  peak,  the 
peak  in  the  power  spectrum  at  the  modulation  frequency,  shows 
no  dip  unlike  any  of  the  optical  spectrum  sidebands.  This  is 
compelling  evidence  that  the  minimization  in  sensitivity  in  the 
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Fig.  11.  Experimentally  measured  (solid  symbols)  and  calculated  (solid 
curves)  values  of  the  modulation  sideband  amplitudes  of  a  normalized  injec¬ 
tion  amplitude  of  £  =  0.06.  This  figure  shows  the  positive  frequency  detuning 
sideband  (+/m ,  solid  circles),  the  negative  detuning  sideband  (— fm ,  solid  tri¬ 
angles),  and  the  corresponding  PI  frequency  (open  squares).  Upper  left  is  for 
the  optical  frequency  component  at  the  frequency  of  the  injected  optical  signal 
(corresponding  to  peak  v\  in  Fig.  5).  Lower  left  is  for  the  component  at  the 
pushed  laser  cavity  resonance  frequency  (corresponding  to  peak  vc  in  Fig.  5), 
one  PI  frequency  below  the  injection  frequency.  Upper  right  is  for  the  side¬ 
band  one  PI  frequency  above  the  injection  frequency  (corresponding  to  peak 
Vi  +  /o  in  Fig.  5),  and  lower  right  is  twice  the  PI  frequency  below  the  injection 
frequency  (corresponding  to  peak  vc  —  /o  in  Fig.  5). 
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Fig.  12.  Mapping  combining  the  calculated  values  of  Fig.  6  and  the  experi¬ 
mentally  measured  values  of  Fig.  8. 
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power  spectrum,  displayed  in  Fig.  10,  arises  from  a  complex 
summation  involving  wave  mixing  in  all  of  the  optical  peaks. 

V.  Discussion  and  Conclusions 

Combining  the  numerical  mapping  results  of  Fig.  5  and  the 
experimental  measurements  of  Fig.  8  yields  the  comparison  map 
of  Fig.  12  that  shows  the  accuracy  of  the  lumped-element  model 
in  describing  the  experimental  results.  Because  it  is  extremely 


difficult  to  determine  the  actual  power  injected  into  the  slave 
laser  through  the  fiber-optic  coupling  path,  this  is  used  as  our 
only  fitting  parameter.  Nevertheless,  it  is  determined  by  fitting 
the  negative  offset  saddle-node  bifurcation  experimental  data  to 
the  model  predictions.  This  then  fully  determines  the  experimen¬ 
tal  map  and  the  agreement  is  excellent.  Further,  the  agreement 
between  data  and  model  in  Fig.  11  shows  that  even  detailed 
features  of  the  system  can  be  fully  recovered  in  quantitative 
agreement.  This  degree  of  accuracy  requires  careful  determi¬ 
nation  of  the  dynamic  parameters  that  go  into  the  model.  The 
parameters  relating  the  real  (gain)  and  imaginary  (refractive  in¬ 
dex)  parts  of  the  complex  gain  to  changes  in  the  carrier  density 
and  circulating  intensity  are  particularly  important  for  setting 
the  positions  of  the  Hopf  and  saddle-node  bifurcations,  which 
strongly  depend  on  the  linewidth  enhancement  factor,  b  [37]. 
The  gain  saturation  factor,  b ',  plays  a  critical  role  in  determining 
the  offsets  of  the  various  minima  in  sensitivity  [31].  The  various 
relaxation  rates  determine  the  magnitude  of  the  PI  frequency 
[38].  Recent  numerical  calculations  show  that  all  of  these  terms 
play  a  role  in  determining  the  positions  of  the  sensitivity  minima 
within  the  overall  mapping  as  a  function  of  injection  amplitude 
and  frequency  offset  [31].  Numerical  work  also  suggests  that 
allowable  values  of  b  for  the  observation  of  the  sensitivity  min¬ 
ima  are  bounded  both  from  below  and  above.  Saturation  effects 
do  not  have  to  be  included  for  the  observation  of  sensitivity 
minima.  Therefore,  they  are  not  determined  by  nonlinear  gain 
effects  such  as  spectral  hole  burning  or  carrier  heating.  It  is  the 
interband  carrier  dynamics,  and  their  effects  on  the  refractive 
index  of  the  gain  medium,  that  principally  determine  when  the 
wave  mixing  of  the  various  optical  frequency  components  will 
lead  to  destructive  interference  and  the  minimization  of  sen¬ 
sitivity  to  perturbations.  Saturation  effects,  and  nonlinear  gain 
effects  in  general,  do  play  a  key  role  in  the  relative  positioning 
of  the  sensitivity  minima  with  respect  to  the  various  bifurcations 
in  the  mapping. 

Because  of  the  excellent  agreement  between  experimental 
data  and  numerical  results,  we  have  confidence  that  the  nu¬ 
merical  simulations  give  an  accurate  representation  of  system 
dynamics,  including  features  that  are  not  easily  accessible  by 
experiment.  Specifically,  the  numerical  simulations  can  be  used 
to  investigate  the  modulation  characteristics  of  the  carrier  den¬ 
sity.  Fig.  13  compares  the  calculated  amplitude  of  the  sidebands 
of  the  PI  frequency  peak  in  the  power  spectrum  with  the  same 
sidebands  in  the  carrier  spectrum.  The  latter  have  been  offset 
by  +17  dB  for  the  negative  sideband  and  by  +18  dB  for  the 
positive  one  from  their  calculated  values  to  show  the  strong 
similarity  between  the  changes  in  the  sideband  amplitudes  in 
the  two  spectra  as  the  offset  frequency  of  the  master  laser  is 
varied.  The  interference  among  the  various  optical  frequency 
components  that  minimizes  the  amplitude  of  the  power  spec¬ 
trum  sidebands  simultaneously  drives  down  the  sensitivity  of 
the  carrier  density  around  these  special  operating  points.  There 
is  minimal  response  of  the  medium  to  the  nonlinear  interac¬ 
tion  of  the  optical  frequency  components  for  frequencies  near 
the  PI  resonance  frequency,  while  the  linear  response  to  the 
induced  current  modulation  is  essentially  unchanged.  Modulat¬ 
ing  the  master  laser  bias  current  simultaneously  changes  the 
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Modulation  Frequency  (GHz) 

Fig.  1 3 .  Calculated  values  of  the  PI  peak  sideband  amplitudes  as  a  function  of 
master  laser  detuning  for  the  intensity  and  carrier  spectra:  (solid  black)  intensity 
upper  sideband,  (solid  blue)  intensity  lower  sideband,  (solid  red)  carrier  upper, 
and  (dashed  white)  carrier  lower.  The  carrier  lower  sideband  has  been  offset  by 
+  17  dB  and  the  carrier  upper  sideband  by  +18  dB  to  emphasize  the  similarity 
with  the  corresponding  intensity  sideband.  Also  shown  is  the  PI  frequency 
(solid  green). 


amplitude  and  frequency  of  the  detuning  so  that  the  operating 
point  would  move  along  a  diagonal  line  in  Fig.  12,  not  a  ver¬ 
tical  line  as  is  the  case  for  temperature  fluctuations.  Because 
the  injected  amplitude  is  changing,  this  perturbation  directly  in¬ 
fluences  the  optical  field  in  (5)  through  a  modulated  injection 
term,  £,  rather  than  through  the  carrier  density  in  (7)  as  is  the 
case  when  the  slave  laser  bias  current  is  modulated.  Past  work 
has  demonstrated  AM-to-FM  conversion  in  the  PI  dynamics 
of  the  optically  injected  semiconductor  laser  [14].  The  slave 
laser  modulation  current  primarily  induces  a  frequency  modu¬ 
lation  and  so  its  influence  on  the  PI  frequency  is  suppressed  at 
these  operating  points.  However,  the  laser  remains  sensitive  to 
fluctuations  in  the  amplitude  of  the  injected  signal. 

These  results  are  interesting  on  two  levels.  On  a  basic  level, 
the  quantitative  agreement  between  theory  and  modeling  allows 
new  tests  of  semiconductor  laser  characteristics  and  dynamic 
parameters  based  on  the  interactions  described  here.  For  in¬ 
stance,  there  is  considerable  current  interest  in  the  dynamics  of 
quantum-dot  lasers  and  the  precision  of  the  techniques  described 
here  will  allow  rigorous  testing  of  proposed  modifications  to  the 
coupled  equation  model  to  characterize  these  novel  structures. 
Further  the  observed  ability  to  reduce  the  sensitivity  of  semi¬ 
conductor  laser  limit-cycle  dynamics  to  system  perturbations 
has  received  very  little  attention  to  date,  and  the  general  charac¬ 
teristics  are  not  well  understood.  Past  work  has  established  that 
the  amplitude  stability  of  free-running  and  stably-locked  semi¬ 
conductor  lasers  can  be  increased  beyond  the  standard  quantum 
limit  [28],  [29].  Here  in  contrast  we  have  demonstrated  that  the 
fundamental  resonance  frequency  in  the  optically  injected  laser 
can  be  made  insensitive  to  fluctuations  in  the  laser  current  and 
operating  temperature  so  that  it  is  only  sensitive  to  the  ampli¬ 
tude  of  the  injected  optical  signal  from  the  master  laser.  This 
leads  to  the  second  interesting  aspect  of  this  paper:  Photonic 
oscillators  have  been  proposed  for  a  variety  of  applications  and 


the  discovery  of  these  operating  points  shows  that  the  optically 
injected  laser  may  be  able  to  display  favorable  phase  noise  and 
jitter  characteristics  to  go  along  with  its  wide  tuning  range  [21], 
[24].  Potentially,  these  operating  points  can  be  complemented 
by  other  nonlinear  dynamics  noise  reduction  techniques  such  as 
those  recently  implemented  in  nanomechanical  oscillators  [39]. 
Now  additional  work  is  needed  to  see  if  these  observations  can 
be  translated  into  improved  performance  of  photonic  devices 
such  as  broadly  tunable,  low-noise,  low-jitter  oscillators. 
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